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1.  Introduction 


The  purpose  of  this  report  is  to  document  the  software  tool  developments  and 
research  findings  made  by  Harry  L.  Hurd  Assoc,  during  the  period  October  i.  1989 
through  May  31,  1992  under  ONR  contract  N00014-S9-C-0310.  This  contract  was  a 
Phase  II  SBIR  award.  This  contract  resulted  from  theoretical  and  simulation  efforts 
performed  under  ONR  contract  N00014-S6-C-0227  and  software  concept  demonstrations 
provided  under  an  ONR  Phase  I  SBIR  award  N00014-SS-C-0683. 

The  primary  technical  topic  addressed  by  this  work  is  that  of  nonstationary 
stochastic  processes.  During  the  past  few  years,  a  deeper  understanding  has  been 
aquircd  concerning  the  nature  of  certain  types  of  nonstationary  processes.  Since  the 
number  of  ways  that  a  process  can  be  nonstationary  are  countless,  particular  classes  of 
nonstationary  processes  have  been  investigated.  Some  of  the  classes  are:  harmonizable. 
periodically  correlated  (cyclostationary).  almost  periodically  correlated,  periodically 
unitary,  almost  periodically  unitary  and  evoliitionarj/.  Specific  e.xamplcs  t)f 
nonstationary  processes  of  interest  to  the  Navy  include  the  radiated  acoustic  signals 
produced  by  cavitating  propellers  and  other  rotating  machinery,  vibration  signals  from 
turbines  or  gearboxes,  meteorological  processes  (exhibiting  daily  or  yearly  nonstationary 
variations)  and  transient  acoustic  signals  produced  by  a  variety  of  sources. 
Nonstationary  processes  arise  in  the  communications  context  where  signals  are  often 
given  a  cyclostationary  structure  in  order  to  facilitate  the  recovery  of  timing 
information  from  the  received  signal.  .Another  class  of  nonstationary  processes  is  given 
by  the  class  of  transient  signals.  It  has  also  been  determined  that  transient  signals  often 
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produce  high  levels  of  spectral  coherence,  and  the  mathematical  reason  for  this 
cjbserva,tion  is  understood. 

This  effort  Wcis  motivated  by  the  recognition  that  the  solution  of  technical 
problems  involving  nonstationary  processes  will  require  a  body  of  theory  and  practical 
tools  that  facilitate  observing  and  characterizing  various  types  of  nonstationary 
processes.  Experimental  tools  in  today's  environment  may  be  taken  to  mean  software 
tools.  Therefore,  this  effort  has  addressed  (1)  continued  research  on  theoretical  issues 
concerning  nonstationary  processes  and  (2)  the  refinement  and  development  (d 
experimental  software  tools  that  will  facilitate  the  investigation  and  characterization  of 
nonstationary  (including  transient)  data. 

The  remainder  of  this  report  contains  (Sec.  2)  a  sketch  of  the  theoretical 
background  that  motivates  our  experimental  methodology,  a  summary  (Sec.  3)  of  new 
theoretical  results  funded  or  partially  funded  by  this  contract,  a  high-level  description 
(Sec.  4)  of  the  workstation  software  package,  called  DSCOH.  that  was  developed  under 
this  contract.  In  addition,  the  use  of  this  software  by  other  Navy  labs  and  projects 
(Sec.  5),  and  the  plan  for  commercialization  (Sec.  6)  will  be  briefly  discussed. 

2.  Theoretical  Background  for  Spectral  Coherence 

Our  general  approach  to  the  development  of  statistical  analysis  tools  and 
algorithms  for  uonstationary  and  transient  processes  is  leased  on  the  theory  of 
harmonizable  stochastic  processes[5].  The  primary  idea  is  that  the  random  amplitudes 
in  the  Fourier  decompositions  of  stationary  processes  are  uncorrelated  or  orthogonal.  In 
contrast,  the  random  amplitudes  in  the  Fourier  decompositions  of  harmonizable 
nonstationary  processes  are  not  generally  orthogonal  in  the  mean-square  sense.  It  is 
f'lear  that  there  may  be  many  ways  in  which  the  random  amplitudes  will  not  be 
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orthogonal  and  hence  many  ways  in  which  a  process  may  be  nonstationary  and  ,et  have 
some  particular  organization  in  the  correlation  of  these  random  amplitudes.  We  now 
discuss  a  little  more  precisely  the  notion  of  harmonizable  processes. 

A  second  order  stochastic  process  X(t)  is  called  harmonizable  [5]  if  it  can  be 
represented  as  a  Fourier  integral  with  respect  to  a  random  measure  Z;  that  is 

X(t)  =  J  exp(iAt)Z(dA)  (1) 

-CO 

where  the  integral  is  in  the  sense  of  quadratic  mean.  Although  there  are  wrakt,!-  notions 
of  harmonizability,  here  we  refer  to  the  strong  notion  for  which  the  spectral  distribution 

r2([a,b)x[c,d))  =  E{[Z(l))-Z(a)][Z(d)-Z(c)]}  (2) 

is  of  bounded  variation  on  RxR.  We  express  this  condition  by 

TO  CO 

J  J  IryldAi.dAo)!  <  x  (3) 

-TO  -00 

and  it  follows  that  the  covariance  function  is  given  by 

R(s,t)  =  E{X(s)X(t)}  =  J  J  exp(isA,-itA2)r7(dAi.dA2).  (4) 

-oc  -OO  "  Lt 

A  second  order  stochastic  process  X(t)  is  said  to  be  (wide-sense)  stationary  when 
the  covariance  function  depends  only  on  the  lag,  i.e.. 

R(s.t)  =  R(s-t)  for  all  s  and  t.  (5) 

Every  mean  square  continuous  stationary  process  has  representation  (1)  where  the 
random  measure  Z  is  orthogonally  scattered  (Z  has  orthogonal  increments)  and  this  will 
occur  if  and  only  if  r2  concentrates  on  the  diagonal.  So  X(t)  being  harmonizable  but 
nonstationary  means  that  Xft)  has  a  decomposition  in  terms  of  complex  exponentials 
exp(iAt)  but  the  increments  of  Z  arc  no  longer  orthogonal.  In  this  sense  harmonizable 
proc('sses  provide  a  natural  framework  for  the  investigation  of  processes  having  spectral 
correlation. 


■A  second  order  stochastic  process  X(t)  is  said  to  be  periodically  correlated  (PC) 


[1],[2],[6]  with  period  T  <  oc  when 


R(s,t)  =  R(s4-T.  t+T)  for  all  s  and  t. 


(6) 


If  X  is  harmonizable.  then  X  is  periodically  correlated  with  period  T  if  and  only 

if  the  support  of  r^  is  contained  in  the  set  S  =  uSj.  where  =  {(Ai,A2);  A^  =  A;  - 

Ic 

2:rk/T},  a  set  of  lines  parallel  to  and  including  the  main  diagonal  as  illustrated  in 
Figure  1.  (See  [G]  for  a  review  of  the  properties  of  harmonizable  PC  processes.) 
Similarly.  X  is  almost  periodically  correlated  if  and  only  if  the  support  of  r^  is 
contained  in  S  =  {(Ai.Ao):  A_,  =  A;  -  ■),  7  6  A  a  countable  set  of  real  numbers  [R2]. 


So  we  .see  that  a  particular  organization  of  the  random  amplitudes  is  that 
associated  with  periodically  correlated  and  almost  periodically  correlated  processes.  To 
paraphrase  the  more  mathematical  statement  above,  a  harmonizable  process  is 
periodically  correlated  if  and  only  if  the  complex  Fourier  coefficients  having  constant 
difference  frequencies  have  non-zero  correlation.  The  pairs  of  frequencies  having 
constant  differences  are  specified  by  a  diagonal  line  in  a  two-dimensional  or  bi-frequency 
plane,  as  illustrated  in  Figure  1.  The  computation  of  diagonal  spectral  coherence  looks 
tor  this  particular  coherence  structure  between  the  complex  Fourier  coefficients 
(amplitudes)  arising  from  discrete  Fourier  transforms  made  on  experimental  time  series. 


2.1  The  Syiectral  Coherence  Computation 

We  shall  now  summarize  the  computation  of  empirical  spectral  coherence  (see 
[Rl])  used  thus  far  in  this  work. 


.A.S  background  for  our  discussion  of  spectral  coherence  we  review  the  notion  of 

coherence  for  complex  Gaussian  random  variables.  Suppose  {Um,ni=l.M}  and 

{\  m-in=l.M}  are  each  sequences  of  independent,  complex  zero  mean  Gaussian  random 
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variables  with  E{|Um|“}='^^'»  E(|Vin|"}=  <^y  ‘ind  E{UmVtn}='Tj^rY.  If  the  values  of  cr^, 
and  (Ty  are  known,  then  a  test  for  >  0  t)e  constructed  from  the  sample  corre¬ 
lation 

t  UmVni-  (~) 


When  (7y  and  ay  are  unknown,  the  normalized  complex  correlation,  called  coherence. 
=  (f\]yl<y\]<yy  may  be  used.  The  usual  estimator  for  7  (actually  l7j~)  is  given  by 


7l" 


M  — ,•) 

I  S  U,nVm|- 

in=  1 _ 

M  --  ,■)  M  i*> 

S  lEtnl-  ^VlVml" 

m  =  1  ni  =  I 


(8) 


X.  R.  Goodman  studied  the  distribution  of  among  other  (piantities  related  to 
complex  Gaussian  random  variables  in  his  di.ssertation  [3].  Goodman  later  suggested  the 
application  of  these  notions  to  testing  for  nonstationarity  [4]  by  interpreting  Um  and 
Vm  to  be  subsets  of  the  discrete  Fourier  transform  vector  {Xj^,k=0.N-l}  where 

%=  i:^X„exp(i2Tkn/N).  (9) 

^  11=0 

.A  test  for  nonstationarity  may  be  expressed  as  a  test  for  a  significantly  large  value  of 
I7I";  or.  in  other  words,  for  correlation  between  the  Um  and  Vm- 

Clearly  there  are  many  ways  to  choo.se  the  {Um}  and  {Vm}  from  {Xj..k=U. 
N-1}.  We  have  applied  this  technique  to  testing  for  the  presence  of  the  PC  property  by 
choosing  {Um}  and  {Vm}  according  to  the  relationship  between  frequencies  permitted 
by  PC  processes.  For  example,  if  {Um}  is  a  contiguous  set  of  elements  of  X.  say  Um  = 
X(j+m)  m=l.  M.  and  we  wish  to  test  for  a  PC  property  at  some  difference  frequency  d. 
then  {V’m}  is  defined  by  Vm  =  X(j-fd-fm)  m=l.  M. 

This  algorithm  can  be  clarified  by  defining  a  two  dimensional  discrete  frequency 
periodogram  from  Xi  .  That  is  we  define 


I 


(10) 


g(N,p.q)  =  X„X 


2tX 


If  the  sequence  Xn  in  (9)  is  PC  with  period  T  (T  an  integer),  then  X,,  is 
liarmonizable  (see  [1])  and  the  support  s<u  ot  the  two-dimensional  spectral  distribution 
is  {(ApA.,)  G  [-t.t):  A.,  =  A^-2Tk/T}.  Further,  the  support  set  Sj.  corresponding  to  A.^ 
Ap27rk/T  transforms  to  2Tq/N  =  2Tp/N-2Tk/T  or  p-q  =  kN/T.  In  other  words,  the  two 
dimensional  discrete  frequency  periodogram  should  be  smoothed  over  lines  of  constant 
difference  frequency  to  search  for  the  presence  of  the  PC  property.  The  spectral 
coherence  calculation  of  Goodman  provides  a  way  of  doing  this  without  knowledge  ot 
the  individual  variances. 


In  order  to  facilitate  this  <l(‘termination  for  an  unknown  time  series  we  compute' 
the  diagonal  spectral  coherence  at  co-ordinate  (p.q)  according  to 


|7(p,q,M)|-  =  - 

’i:  |X 


M-l  -c - 

v  V  Y  ( 

m=()  P  + 


(11) 


We  compute  these  values  in  a  sepiare  array  (over  p.q)  and  plot  the  values  that  exceed  a 
threshold.  -A.  useful  threshold  computation  is  that  given  by  Goodman  for  the  null  distri¬ 
bution  for  It!"  under  the  assumption  that  {Cm}  <^Ild  {Vpi}  are  complex  Gaussian,  arise 
from  disjoint  frequency  sets  and  the  variances  E{|Uni|”}  =  E{|V,t,|“}  =  ay  are 

constant  with  respect  to  m.  In  this  case 

trllTt'  >  Iron  = 


.3.  N<'w  Theoretical  Results 

Theoretical  research  on  non-stationary  processes  is  yielding  fundamentally  new 
information  about  the  structure  of  perioflically  and  almost  periodically  correlated 
processes;  specifically,  we  have  shown  how  groups  of  unitary  operators  are  intimatc'ly 
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connected  with  these  processes.  These  research  findings  will  permit  the  solution  of 
[jroblems  that  were  difficult  when  working  only  with  the  correlation  theory.  In  addition, 
these  findings  provide  representations  that  are  extendable  to  multiple  dimensions,  that 
is.  to  random  fields.  These  findings  will  also  clarify  the  connection  between  periodically 
correlated  processes  and  wavelet  representations  of  processes.  The  research  documented 
in  the  following  manuscripts  has  been  partly  or  entirely  supported  by  this  contract. 

[Rl]  Hurd.  H.L.,  and  Gerr,  N..  Graphical  Methods  for  Determining  the  Presence  of 
Periodic  Correlation.  J.  Time  Series  .Anal..  12,  No.  4.  337-350.  1991. 

[R2]  Hurd.  H.L..  Correlation  Theory  of  Almost  Periodically  Correlated  Processes,  .J. 
Multivariate  .4nal..  37,  No.  1,  24-45,  1991. 

[R3]  Hurd.  H.L..  The  Unitary  Operator  of  Periodically  Correlated  Processes.  Invited 
paper  given  at  the  Hampton  Univ.  Conference  on  Nonstationary  Stochastic 
Processes.  .4ug.  1  and  2.  1991. 

[R4]  Hurd,  H.L..  Almost  Periodically  Unitary  Stochastic  Processes.  Tech.  Rept.  No. 
284,  Center  for  Stochastic  Processes.  Dept,  of  Statistics.  University  of  NC  at 
Chapel  Hill,  February  1990.  To  appear  in  Stochastic  Processes  and  their 
.Application. 

[R5]  Houdre',  C.H..  On  the  Spectral  Strong  Law  of  Large  Numbers  and  the  Pointwise 
Ergodic  Theorem  in  L'-'’ .  To  appear  in  .Annals  of  Probability. 

[R6]  Hurd.  H.L.  and  Leskow.  .1..  Strongly  Consistent  and  Aysmptotically  Norm.al 
Estimation  of  the  Covariance  for  Almost  Periodically  Correlated  Processes.  Tech. 
Rept.  No.  308.  Center  for  Stochastic  Processes.  Dept,  of  Statistics.  University  of 
NC  at  Chapel  Hill.  Sept.  1990.  To  appear  in  Statistics  and  Decisions. 

[R7]  Cambanis.  S.,  Houdre.  C.H..  Hurd,  H.L.,  and  Leskow.  J..  Laws  of  Large  Numbers 
for  Periodically  and  Almost  Periodically  Correlated  Processes.  Tech.  Rept.  No. 
334.  Center  for  Stochastic  Processes.  Dept,  of  Statistics.  University  of  NC  at 
Chapel  Hill.  March  1991.  Submitted  for  publication. 

iR8]  Hurd.  H.L.  and  Leskow.  .1.,  Estimation  of  the  Fourier  Coefficient  Functions  and 
their  Spectral  Densities  for  o-Mixing  Almost  Periodically  Correlated  Processes. 
Tech.  Rept.  No.  330,  Center  for  Stochastic  Processes.  Dept,  of  Statistics. 
University  of  NC  at  Chapel  Hill,  March  1991.  To  appear  in  Statistics  and 
Probability  Letters. 

[R9]  Houdre',  C.H..  Some  Recent  Results  on  the  Sampling  Theorem..  To  appear  in  1992. 

[RIO]  Hurd.  H.L.  and  Mandrekar.  \’..  Spectral  Theory  of  Periodically  and  Quasi- 
Periodically  Stationary  SoS  Sequencess.  Tech.  Rept.  No.  349.  Center  for 
Stochastic  Processes,  Dept,  of  Statistics.  UNC  at  Chapel  Hill.  Sept.  1991. 
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[Rll]  Hurd,  H.L.  and  Russek.  A..  Stepanov  Almost  Periodically  Correlated  and  Almost 
Periodically  Unitary  Processes.  Tech.  Rept.  No.  368,  Center  for  Stochastic 
Processes,  Dept,  of  Statistics.  UNC  at  Chapel  Hill.  August.  1992. 

[R12]  Hurd,  H.L.  and  Kallianpur.  G..  Periodically  Correlated  and  Periodically  Unitary 
Processes  and  their  Relationship  to  ’L.yfO.T] -  Valued  Stationary  Sequences.  In 
Nonstationary  Stochastic  Processes  and  Their  Application,  .I.C. Hardin  and 
A.G.jMiamee  Ed..  World  Scientific  Publishing  Co..  1992. 

[R13]  Houdre,  C.H.,  On  the  Irregular  Sampling  Theorem.  In  progress. 


[R14]  Hurd.  H.L.  and  Russek.  A..  Almost  Periodically  Correlated  Processes  on  LCA 
Groups,  Tech.  Rept.  No.  369.  Center  for  Stochastic  Processes.  Dept,  of  Statistics. 
UNC  at  Chapel  Hill.  1992. 

[R15]  Bloomfield.  P..  Hurd.  H.L..  and  Lund.  R..  Periodic  Correlation  in  Meteorlogical 
Time  Series.  Presented  at  the  3th  International  Meeting  on  Statistical 
Climatology,  Toronto.  Quebec.  .June  22-26.  1992. 

[R16]  Bloomfield.  P.,  Hurd.  H.L..  and  Lund.  R..  Periodic  Correlation  in  Stratospheric 
Ozone,  Tech.  Rept.  No.  2076.  Institute  for  Statisticss.  Dept,  of  Statistics.  L^NC  at 
Chapel  Hill.  June  15.1992. 


In  addition,  the  following  application  oriented  reports  have  been  produced. 


[Al]  Lund.  R..  A  Sensitivity  Study  of  Goodman's  Coherence  Statistic  Via  Simulations. 
Technical  Report  by  Harry  L.  Hurd  Associates,  on  Contract  N00014-S9-C-0310  . 
September  15.  1991. 

[.42]  Lund,  R..  The  Kernel  Method  of  Density  Estimation:  A  Summary,  Technical 
Report  bv  Harrv  L.  Hurd  .A.ssociate.s.  on  Contract  N00014-S9-C-0310  .  October  14. 
1991. 

[.\3|  Jones.  C.,  Spectral  Coherence  of  Harmonics.  Consultant’s  report  submitted  to 
Harry  L.  Hurd  .Associates,  on  Contract  N00014-S9-C-0310  ,  November  5,  1991. 

[.44]  Jones,  C..  PC  Signal  Models.  Consultant's  report  submitted  to  Harry  L.  Hurd 
.Associates,  on  Contract  N0()014-S9-C-0310  .  December  9.  1991. 

[A5]  Jones.  C..  Helicopter  Data  Analysis.  Consultant's  report  submitted  to  Harry  L. 
Hurd  Associates,  on  Contract  N00014-S9-C-0310  ,  February  3,  1992. 

[.46]  Jones.  C..  DSCOH  User's  Ma.nual.  Technical  Report  bv  Harrv  L.  Hurd  .Associates, 
on  Contract  N00014-S9-C-0310  .  July  31.  1992. 

[.47]  Hurd.  H.  and  Jones.  C..  DSCOH  Technical  Reference  Manual.  Technical  Report 
by  Harry  L.  Hurd  .Associates,  on  Contract  N00014-S9-C-0310  .  July  31.  1992. 
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[AS]  Greineder.  S.  and  Hurd.  H.,  Investigation  of  Diagonal  Spectral  Coherence  as  a 
Class  Discriminant  Function.  XUSC  Technical  Report(SECRET),  in  preparation. 


4.  Workstation  Software 

The  dev’elopmental  part  of  the  effort  on  contract  X00014-S9-C-0310  has  focused 
on  the  development  of  workstation  .software  to  facilitate  the  examination  and 
characterization  of  nonstationary  data  using  spectral  coherence.  We  have  continued  to 
improve  the  software  functionality  and  u.ser  interface  over  the  span  of  this  project. 
Early  in  the  project,  each  basic  function  was  implemented  as  a  stand-alone  program. 
.\t  the  next  stage,  the  separate  inograms  were  integrated  into  one  master  progrant 
called  TD.  Then  versions  of  the  program,  called  TDP.  were  built  that  could  use  array 
processor  cards  in  both  the  IDM-PC  family  or  in  the  MacII  family.  In  the  latter  part  ol 
the  project,  we  have  concentrated  on  the  IBM-PC  family  because  of  programming  ease 
and  because  the  S04S6  chip  gives  ample  speed  for  the  basic  functions,  even  without  the 
array  processor.  The  latest  version,  called  DSCOH.  contains  a  greatly  improved  user 
interface,  and  is  considered  the  final  i)roduct  of  the  contract. 

We  shall  first  summarize  the  basic  functionality  of  DSCOH.  followed  by 
functionality  provided  by  other  programs  developed  under  this  contract.  These  other 
programs  include  the  earlier  TD.  TDP  and  TDNEW  programs  or  by  more  recent  stand¬ 
alone  programs. 

4.1  DSCOH  Functionality 

(a)  Plotting  of  time;  series. 

(b)  Periodograms.  Periodograms  are  plots  of  the  squared  magnotude  ol  sample 

Fourier  transforms. 

(c)  Spectrograms,  Spectrograms  are  waterfall  or  gram  displays  made  from 

successive  periodograms  using  various  amplitude  encoding  and  scaling  methods. 
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(il)  Spectral  colK'n-iua*  iiuaee  (lisf)iavs.  The  SC  images  are  cpniputecl  aeeordiiit>;  ti, 
Eq.(ll)  and  grev-srale  encoded.  This  display  eoinhiiies  five  visual  outputs  bas('d  on 
the  FFT  of  a  segnu'iit  of  data.  These  include  th«'  SC  iinag(>,  a  time  serit's  plot  ot 
the  current  data  sement.  plots  of  the  usual  periodogram  and  the  phase  spectrum, 
and  two  auxilliarv  plots.  The  atixilliary  plots  may  lx-  chosen  in  pairs  from  (dther  of 
'he  following  two  sets. 

(i)  Plots  supportinu  tin*  detection  ot  periodic  coirelation. 

These  plots  both  are  ba.sed  on  tin*  fact  that  the  support  set  for  a  PC  proci'ss  is  a 
collection  of  straight  lines  jiarallel  to  the  main  iliagonal.  The  fn'st  plot  in  this 
pair  may  be-  motivated  by  the  thought  experiment  of  increasing  M  to  be  as 
large  as  permitted  so  that  only  one  value  ot  |-(p.q.M)|”  is  determined  for  each 
separation  d  from  the  main  diagonal.  In  ti-rms  of  |'(p-q.Ml|“  compute 
|-(0.d.X  '1)1~  where  X  is  the  diiiKUision  of  the  Fourier  transform  vc'ctor.  The 
implementation  is  btised  on  an  e<puvalent  but  simpler  computation.  I'onsistinu 
of  the  DFT  of  the  s([uare.s  of  the  input  time  series.  The  second  jilot  is  based  on 
the  incoherent  averaging  of  all  the  values  of  i-;.(p.(i.M)|“  from  a  cohertuice 
surface  for  a  fixed  t-alne  of  d.  Symbolically,  wt'  computt'  the  incoherent  average 
according  to 

^(d.XI)  =  JL  f  hdpM.  p.M+a.  Mir'  (12) 

^  +  ^  p=() 

where  L  =  [(X’  — l-d)/Mj  and  IxJ  denotes  the  gnsitest  integer  in  x. 

(ii)  Plots  siirijioi'tiii'j.'  the  detection  qt  transients. 

The  first  plot  in  this  pair  is  a  plot  ot  the  TR-1  statistic,  which  is  determim'd  bv 
the  h’vel  of  coherence  obtained  by  traversing  the  spectral  coherence  images  on  ;i 
diagonal  path  of  difference  frequency  d.  To  be  more  precise,  the  TR-1  ordinates 
are  given  by  |-;dk.k+d.M)|  as  a  function  of  k.  It  is  thought  that  small  values  of 
d  are  appropriate.  The  second  plot  in  this  pair  is  a  plot  of  the  TR-2  statistic, 
which  is  determined  })y  the  ratio  ot  a  weighted  coherent  sum  of  the  Fourier 
coefficients  to  the  incoherent  sum 

X-1  ~  .,  X-1  ~ 

=  \T  X.  ,  ,  exp(iAj)|-/  V  L.r  (13) 

j=l)  j=0 

The  algorithm  is  motivated  by  the  ob.servation  that  if  the  siiectrum  of  a 
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transient  signal  has  a  linear  phase  function  exp(-iAj)  in  the  neighl)orhood  of  the 
maximum  spectral  levels,  then  the  statistic  ;/(k)  will  be  large  (>1)  when  A  is 
chosen  to  make  the  terms  in  the  sum  add  coherently.  A  disadvantage  of  this 
statistic  is  that  multiple  values  of  A  may  have  to  be  tried.  In  contrast,  the 
algorithm  TR-1  is  much  more  robust  because  its  output  is  independent  of  A  for 
a  linear  phase  function  o(j)=Aj. 

(e)  Correlation  and  cepstral  ]>lors. 

■Autocorrelation  and  cepstral  plots  are  provided  to  assist  in  det(umining  the 
presence  of  multi-path  signals. 

4.2  Other  Functionality 

(a)  PC-urarns.  .4  PC-gram  is  formed  in  the  same  manner  as  the  spectrogram  but 
from  sticcessive  traces  from  ou<‘  of  the  auxilliary  plots  supporting  detection  of  PC 
processes.  Either  auxilliary  plot  may  Ite  chosen  as  the  source  of  the  PC-gram.  The 
region  of  the  SC  image  contributing  to  the  incoherent  sum  maj"  be  restricted  by  the 
user. 

(b)  T-grams.  T-grams  are  grams  made  from  successive  scans  of  the  TR-1  transient 
detection  statistic. 

(c)  Bandshifting.  filtering  and  decimation.  This  basic  function  permits  the 
investigator  to  produce  data  files  containing  filtered,  bandshifted  (to  center- 
frequency  of  zero)  and  decimated  sub-bands  of  time  series  from  a  specified  input 
file. 


•5.  L  se  of  Workstation  Softwan*  bu  Navv  Labs  and  Projects 


In  the  course  of  this  project,  c'xperimental  versions  of  the  software  have  been 
installed  at  various  Xavy  labs.  This  provided  X’avy  scientists  with  capabilities  for 
analyzing  nonstationary  signals  using  spectral  coherence,  and  also  provided  feedback 
from  which  improvements  were  based.  In  particular,  versions  of  the  software  have  been 
provided  to  NUWC  (formerly  NUSC).  XCCOSC  XRaD  (formerly  XOSC)  XSWC  and 
DTRC.  These  interactions  have  lead  to  contractural  support  for  the  investigation  of 
spectral  coherence  methods  in  specific  aiiplications.  Support  for  this  application  work  is 
coming  from  XUWC.  XADC  and  XCCOSC  (the  DARPA  DAXTES  program).  Final 
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versions  of  the  software  have  been  deliven'd  to  XOSC  and  NUSC  as  called  out  in  the 
contract. 


6.  Commercialixation  of  the  Workstation  Software 


The  DSCOil  software  package  will  Ije  offered  for  public  sale  in  the  near  future. 
The  capabilities  offered  b}’  DSCOH  will  permit  an  investigator  to  perform  basic 
analyses  of  signals  that  are  thought  to  have  diagonal  spectral  coherence.  Further  details 
on  the  DSCOH  program  may  be  found  in  the  DSCOH  User's  Manual  and  the  DSCOH 
Technical  Reference.  Manual. 
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Figure  1— Spectral  Mass  Locations  for 
Ilarmonizable  Periodically  Correlated  Processes 
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